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Abstract 

We study the hydrogen and muonic hydrogen within an effective fieid theory frame- 
work. We perform the matching between heavy baryon effective theory coupled to pho- 
tons and leptons and the relevant effective field theory at atomic scales. This matching 
can be performed in a perturbative expansion in a, l/m p and the chiral counting. We 
then compute the 0(mf.a 5 /mH x logarithms) contribution (including the leading chiral 
logarithms) to the Hyperfine splitting and compare with experiment. They can explain 
about 2/3 of the difference between experiment and the pure QED prediction when 
setting the renormalization scale at the p mass. We give an estimate of the matching 
coefficient of the spin-dependent proton-lepton operator in heavy baryon effective the- 
ory. 
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1 Introduction 



Years have passed since the advent of QCD. After numerous attempts to understand QCD 
by using several models, more studies now move towards trying to parameterize the QCD 
properties in a model independent way with the help of different systematics that are usu- 
ally highlighted by the specific kinematic situation under study. One could hope that this 
approach may bring some light in the understanding of QCD or a least to provide some 
consistency check between different models. Therefore, it becomes important to be able to 
relate as many observables as possible in a model independent framework. Effective field 
theories (EFT's) may play an important role in this approach. 

Within the above philosophy, the study of hydrogen (ep) and muonic hydrogen ((J,p), 
in particular of the high precision measurement of different splittings, can provide accurate 
determinations of some hadronic parameters related with the proton elastic and inelastic 
electromagnetic form factor, like the proton radius and magnetic moment, polarization ef- 
fects, etc.... 

In the ep and up we are basically testing the proton with different probes (e, /i, 7). 
They correspond to the simplest possible probes since they are point-like particles and the 
interaction is perturbative (the analogy with deep inelastic scattering is evident and it has 
already been used since long ago [1, 2, 3, 4] in order to obtain some of these hadronic 
parameters from dispersion relations). They also provide the first natural step towards more 
complicated systems like exotic or heavy (muonic) atoms. 

The ep and \xp systems are, in a first approximation, states weakly bound by the Coulomb 
interaction and their typical binding energy and relative momentum are, E ~ m e ^a 2 and 
|p| ~ m e ^a, respectively. We will switch off the weak interactions in this work. Therefore, 
the ep and [ip systems become stable and C, P and T are exact symmetries of these systems. 
In any case, several different scales are involved in their dynamics: 

For the ep system they are: . . ., m e a 2 , m e a, m e , Am = m n — m p , A = uia — m p , m n , 
m p , m p , A x , . . . that we will group and name in the following way: 

• m e a 2 : ultrasoft (US) scale. 

• m e a: soft scale. 

• Me» = me " lp Am = m n — m v , m e : hard scale. 

• m M , A = m& — m p , m n : pion scale. 

• m p , m p , A x : chiral scale. 

For the \ip system they are: . . ., m^a 2 , m^a, m M , Am = m n — m p , m e , A = m& — m p , 
m n , m p , m p , A x , . . . that we will group and name in the following way: 

• m^a 2 : US scale. 

• Am = m n — m p , m e , m^a: soft scale. 

• /Vp = m^+m ' "V> A = m A - m p , m w : hard/pion scale. 
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• m p , m p , A x : chiral scale. 

By doing ratios with the different scales, several small expansion parameters can be built. 
Basically, this will mean that the observables, the spectrum in our case, can be written, up 
to large logarithms, as an expansion, in the case of the ep, in a, ^ and and in the 
case of the /ip, in a and ^f-. It will also prove convenient sometimes to use the reduced 
mass /i M ( e )p, since it will allow to keep (some of) the exact mass dependence at each order in 
a. In order to be more precise, the ep energy will be expanded in the following way (up to 
logarithms): 

E(ep) = -^(l + c 2 a 2 + c 3 a* + •••), (1) 

where 

(2) 

and cfrfi are functions of dimensionless quantities of 0(1) like jj, ep /m e , m^/m^, etc. 
For the \xp things work analogously: 

2 

E(np) = -^-(1 + c ± a + c 2 a 2 + c 3 a 3 + •••), (3) 
where c\ does not depend on hadronic quantities, only on m^a/m e , and (n > 2) 

where are functions of dimensionless quantities of 0(1) like n w /m^ m^/m^, m p a/m e , 
etc. 

Let us stress that the coefficients c n can be expanded in the ratio i.e. in the 
chiral/heavy-baryon expansion (m p should also be understood as A x ). 

In order to disentangle all the different scales mentioned above it is convenient to use 
EFT's. In order to obtain the relevant one for these systems we first need to decide what 
are the degrees of freedom we want to describe. In our case we want to describe the ep 
and up systems at ultrasoft or smaller energies. Therefore, degrees of freedom with higher 
energies can (and will) be integrated out in order to obtain the EFT to describe these 
systems. One EFT that fulfills this requirement is potential NRQED (pNRQED) [5, 6] (for 
some applications see [7] and see also [8, 9]). pNRQED appears after integrating out the 
soft scale from NRQED [10] and it shares some similarities with the approach followed in 
Ref. [11]. We will obtain pNRQED by passing through different intermediate effective field 
theories after integrating out different degrees of freedom. The path that we will take is the 
following (in some cases, instead of this chain of EFT's one can use dispersion relations, or 
direct experimental data, in order to obtain the matching coefficients): 

HBET -> (QED) -> NRQED -> pNRQED. 

This way of working opens the possibility to compute the observables of atomic physics with 
the parameters obtained from heavy baryon effective theory (HBET), which is much close to 
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QCD since it incorporates its symmetries automatically, in particular the chiral symmetry. 
Besides, it is the matching with HBET that will allow to relate the matching coefficients used 
for ep with the ones used in up. HBET [12] describes systems with one heavy baryon: the 
proton, the neutron or also the delta [13] at the pion mass scale. The chiral scale explicitly 
appears in the Lagrangian as an expansion in 1/ A x , l/m p and any other smaller scale remains 
dynamical in this effective theory. In short, HBET is a EFT defined with an UV cut-off v 
such that i/< A x but larger than any other dynamical scale in the problem. 

In the up, NRQED appears after integrating out the hard scale, whereas in the ep, 
NRQED appears after integrating out the pion and hard scales. In this last case one could 
pass through an intermediate theory (QED) defined by integrating out the pion scale and 
profit from the fact that pion and hard scales are widely separated. Nevertheless, we will do 
the matching here in one step for simplicity. 

pNRQED is obtained after integrating out the soft scale. We refer to [5, 6] for further 
details. 

The above methodology allows to compute (or parameterize in a model independent way) 
the coefficients c's in a systematic expansion in the different small parameters on which these 
systems depend. 

It is the aim of this paper to use this procedure for the computation of the leading- 
logarithm hadronic contributions to the hyperfine splitting for both the ep and fip. This 
means to compute the spin-dependent piece of C3 with 0((m e /m p ) 2 x logarithms) and 
0((m M /m p ) 2 x logarithms) accuracy for the ep and \xp respectively. 

2 Effective Field Theories 

In this section, we will consider the different EFTs that will be necessary for our calculation. 
2.1 HBEFT 

Our starting point is the SU(2) version of HBEFT coupled to leptons, where the delta is 
kept as an explicit degree of freedom. The degrees of freedom of this theory are the proton, 
neutron and delta, for which the non-relativistic approximation can be taken, and pions, 
leptons (muons and electrons) and photons, which will be taken relativistic. 

Our first aim will be to present the effective Lagrangian of this theory. It corresponds 
to a hard cut-off jj, « m p , A x and much larger than any other scale in the problem. The 
Lagrangian can be split in several sectors. Most of them have already been extensively 
studied in the literature but some will be new. Moreover, the fact that some particles will 
only enter through loops, since only some specific final states are wanted, will simplify the 
problem. The Lagrangian can be structured as follows 

£-HBET = £7 + £-1 + £tt + £-Itt + £(N,A) + £{N,A)l + £(N,A)k + £{N,A)lin (5) 

representing the different sectors of the theory. In particular, the A stands for the spin 3/2 
baryon multiplet (we also use A = — m p , the specific meaning in each case should be 
clear from the context). 



3 



The Lagrangian can be written as an expansion in e and l/m p . Our aim is to obtain 
the hyperfine splitting with 0(mf.a: 5 /m 2 , x (mm 9 ,m A, lnmjj) accuracy, where m q stands 
for the mass of the light, u and d (or s), quarks and for the mass of the lepton (the 
leading order contribution to the hyperfine splitting reads Ep = (8/3)c^mf.a 4 /m p , where 

Cp^ is defined in Eq. (10)). Therefore, we need, in principle, the Lagrangian with 0(l/m 2 ) 
accuracy. Let us consider the different pieces of the Lagrangian more in detail. 



The photonic Lagrangian reads (the first corrections to this term scale like a 2 /m 



C, = -\F^F, U . (6) 

The leptonic sector reads (iD^ = id^ — eA^) 

C^^UtiP-m^k, (7) 



where i — e, fj,. We do not include the term 

egi 



m p 



since the coefficient is suppressed by powers of a and the mass of the lepton. Therefore, it 
would give contributions beyond the accuracy we aim. In any case, any eventual contribution 
would be absorbed in a low energy constant. 

The pionic Lagrangian C n is usually organized in the chiral counting. From the analysis 
of the next section we will see that the free pion propagators provide with the necessary 
precision. Therefore, we only need the free-particle pionic Lagrangian: 

= (0„7r + )(0"7r-) - mln+n~ + i(^7r°)(^7r°) - ^mfrV . (9) 

The one-baryon Lagrangian £(jv,a)tt is needed at 0(l/m 2 ). Nevertheless a closer inspec- 
tion simplifies the problem. A chiral loop produces a factor 1/(AttF ) 2 ~ 1/ra 2 . Therefore, 
the pion-baryon interactions are only needed at 0(m n ), the leading order, which is known 
[12, 13]. 1 For the explicit expressions we refer to these references. 

Therefore, we only need the one-baryon Lagrangian C(n,a) at 0(l/m 2 ) coupled to elec- 
tromagnetism. This would be a NRQED-like Lagrangian for the proton, neutron (of spin 
1/2) and the delta (of spin 3/2). The neutron is actually not needed at this stage. The 
relevant term for the proton reads 

f d 2 D 4 c (p) c (p) 1 

«W> = K\ + ^ + g£ - ^ B - ieZ p ^ a . (D p x E - E x D p ) | JV, , 

(10) 



1 Actually terms that go into the physical mass of the proton and into the physical value of the anomalous 
magnetic moment of the proton n p = Cp^ — 1 should also be included (at least in the pure QED computations) 
and it will be assumed in what follows. For our computation these effects would be formally subleading. In 
any case, their role is just to bring the bare values of mo and no to their physical values. Therefore, once 
the values of m p and /i p are measured by different experiments, they can be distinguished from the effects 
we are considering in this paper. 
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where iD® = id + Z p eA°, iD p = i'V — Z p eA. For the proton Z p = 1. We have not included 
a term like 

(p) 

%Vt[V-E]iV (11) 
ritp 

We could have done so but it may also be eliminated by some field redefinitions. In any case 
it would give contribution to the spin-independent terms so we will not consider it further 
in this work. 

As for the delta (of spin 3/2), it mixes with the nucleons at 0(l/m p ) (0(l/rrip) are not 
needed in our case). The only relevant interaction in our case is the p-A + -7 term, which is 
encoded in the second term of 

5£ {NA) = T^ido - A)T + ^ (lV<$> • B rf/^N + h.c) , (12) 

where T stands for the delta 3/2 isospin multiplet, N for the nucleon 1/2 isospin multiplet 
and the transition spin/isospin matrix elements fulfill (see [14]) 

4m = - -^ fc )> ^tm^rn = - ( 13 ) 

The baryon-lepton Lagrangian provides new terms that are not usually considered in 
HBET. The relevant term in our case is the interaction between the leptons and the nucleons 
(actually only the proton): 

S£(n,a)i = \ E $%n p1 °n p laoh + \ E ct; R N P i 3 isN p lamk ■ (14) 

m p i m p i 

The above matching coefficients fulfill c^ R = R and c^ l R = R up to terms suppressed by 
mijirip, which will be sufficient for our purposes. 

Let us note that with the conventions above, N p is the field of the proton (understood 
as a particle) with positive charge if li represents the leptons (understood as particles) with 
negative charge. 

This finishes all the needed terms for this paper, since the other sectors of the Lagrangian 
would give subleading contributions. 



2.2 NRQED(^) 

In the muon-proton sector, by integrating out the m n scale, an effective field theory for 
muons, protons and photons appears. In principle, we should also consider neutrons but 
they play no role at the precision we aim. The effective theory corresponds to a hard cut-off 
v « and therefore pions and deltas have been integrated out. The Lagrangian is equal 
to the previous case but without pions and deltas and with the following modifications: 
C\ — > C e + C^ R ^ and C(n,a)i ~^ ^Ne + ^n^\ where it is made explicit that the the muon 
has become no relativistic. Any further difference goes into the matching coefficients, in 
particular, into the matching coefficients of the baryon-lepton operators. In summary, the 
Lagrangian reads 

£-NRQED{n) = Cry + C e + C^ R ^ + Cn + £jVe + C^N^ i (^) 
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where 

f D 2 D 4 r (p) r M 1 

Cr - + + + B+fcZ,^, ■ (D, x E — E x D ( ,)j,„ (16) 

and 

£j™ = C ^fKN p 1% - 4&N}aN, , (17) 

with the following definitions: iD® = id — Z^eA , iD^ = + Z^eA and Z^ = 1. C e stands 
for the relativistic leptonic Lagrangian (7) and C Ne for Eq. (14), both for the electron case 
only. 

A term of the type 

' !Jle l e a^l e F^ (18) 



is not taken into account because of the same reason as in Sec. 2.1. 



2.3 QED(e) 

After integrating out scales of 0(m n ) in the electron-proton sector, an effective field theory 
for electrons coupled to protons (and photons) appears. Again, we should also consider 
neutrons but they play no role at the precision we aim. This effective theory has a cut-off 
v « and pions, deltas and muons have been integrated out. The Lagrangian reads 

£-QED(e) = C y + C e + Cn + ^TVe • (19) 

This Lagrangian is similar to the previous subsection but without the muon. 



2.4 NRQED(e) 

After integrating out scales of 0{m e ) in the electron-proton sector, we still have an effective 
field theory for electrons coupled to protons and photons. Nevertheless, now, the electrons 
are non-relativistic. The Lagrangian is quite similar to the one in subsec. 2.2 but without a 
light fermion and with the replacement fi — > e. The Lagrangian reads 

C N RQED { e) = + Cf^ + C N + ■ (20) 



2.5 pNRQED 

After integrating out scales of 0(mi.a) one ends up in a Schrodinger-like formulation of 
the bound-state problem. We refer to [5, 6] for details. The pNRQED Lagrangian for the 
ep (the non-equal mass case) can be found in Appendix B of the second Ref. in [6] up 
to 0(ma b ). The pNRQED Lagrangian for the up is similar except for the fact that light 
fermion (electron) effects have to be taken into account. The explicit Lagrangian and a more 
detailed analysis of this case will be presented elsewhere. For the purposes of this paper, we 
only have to consider the spin-dependent delta potential: 

5 y = 2 ^NR S 2 6 (3) {r y (21) 

which will contribute to the hyperfine splitting. 
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3 Form Factors: definitions 

It will turn out convenient to introduce some notation before performing the matching be- 
tween HBET and NRQED. We first define the form factors, which we will understand as 
pure hadronic quantities, i.e. without electromagnetic corrections. 

Our notation is based on the one of Ref. [15]. We define J M = J2i QiQil^Qi where % — u,d 
(or s). The form factors are then defined by the following equation: 



{p' , s\J»\p, s) = u(p') 



Fi{q)Y + iF 2 (q 



u(p) , (22) 



2m 

where q = p' — p and iq, F 2 are the Dirac and Pauli form factors, respectively. The states 
are normalized in the following (standard relativistic) way: 

(p', X'\p, A) = (27r) s 2p°* S (p' - p)5y X , (23) 

and 

u(p, s)u(p, s) = ^ — 2 — ' ^ > 

where s is an arbitrary spin four vector obeying s 2 = — 1 and P ■ s = 0. 
The form factors could be (analytically) expanded as follows 

F l {q 2 )=F l + ^- 2 F' l +... (25) 
m z 

for very low momentum. Nevertheless, we will be interested instead in their non-analytic 
behavior in q since it is the one that will produce the logarithms. 
We also introduce the Sachs form factors: 

G E (q 2 ) = F 1 (q 2 ) + ^F 2 (q 2 ), G M (q 2 ) = iq(g 2 ) + F 2 {q 2 ). (26) 
We will also need the forward virtual-photon Compton tensor 

= iJd A xe i(1 - x (p,s\TJV{x)r{0)\p,s) , (27) 
which has the following structure (p — q • p/m): 

T^= (-9 >lu + qJ f)S 1 (p,q 2 ) 

+ ^ p (p"-^)(p u -^q l ')S 2 (p,q 2 ) 

-±- v ^q p s a A x {p,q 2 ) 
-^3 e^q p {{m p p)s a - (q ■ s) Pf7 )A 2 (p, q 2 ), (28) 

depending on four scalar functions. It is usual to consider the Born approximation of these 
functions. They read 

QBorn/ 2x _ 4mgq 2 F?(q 2 )-(q 2 ) 2 Fj(q 2 ) , . 

^2 {PiQ ) ~ Z (2m p p)M<? 2 ) 2 ' ^ ) 



In the following section, we will use the results of Ji and Osborne [16]. Their notation 
relates to ours in the following way (for the spin-dependent terms): S{° = Ai/m 2 and 
5*2 — A 2 /m p i . Note however that Af orn above is different from the Si definition in Ref. [16] 
by the F% term. 

4 Matching 

The matching between HBET and NRQED can be performed in a generic expansion in 
l/m p , l/m M and a. We have two sort of loops: chiral and electromagnetic. The former 
are always associated to 1/(47tFq) 2 factors, whereas the latter are always suppressed by a 
factors. Any scale left to get the dimensions right scales with m n . In our case we are only 
concerned in obtaining the matching coefficients of the lepton-baryon operators of NRQCD 
with 0(a 2 x (lnm 9 , In A, lnm/J) accuracy. Therefore, the piece of the Lagrangian we are 
interested in reads 

SC = ^fNlN p l\k - ^Nl<rN p ijah , (33) 

where we have not specified the lepton (either the electron or the muon). In what follows, we 
will assume that we are doing the matching to NRQED(/x). Therefore, we have to keep the 
whole dependence on mijm^. The NRQED (e) case can then be derived by expanding m e 
versus m^. In principle, a more systematic procedure would mean to go through QED(e). 
Nevertheless, for this paper, as we do it will turn out to be the easier. 

In principle, the contributions scaling with are the more important ones. Never- 

theless, they go beyond the aim of this paper. This is specially so as far as we are only 
interested in logarithms and the spin-dependent term c^* NR . Its general expression at 0(a 2 ) 
reads (an infrared cutoff larger than mia is understood and the expression for the integrand 
should be generalized for an eventual full computation in D dimensions) 

*™ - -t / S^h^m l AiiU ^ + 2k2) + ' ik2 ^t Mk< " k2) ) • (34) 

consistent with the expressions obtained long ago like in Ref. [3]. This expression has 
been obtained in the Feynman gauge. It correctly incorporates the whole dependence on 
the lepton mass. Therefore, the same expressions are valid for the hydrogen and muonic 
hydrogen. 

In principle one should also consider contributions with one, two or three electromagnetic 
current insertions in the hadronic matrix elements instead of only two as in Eqs. (27) and 
(34). Nevertheless, only the above contribute to the order of interest. 

Within the EFT framework the contribution from energies of 0(m p ) or higher in Eq. (34) 
are encoded in c^ l R ~ c^ R (analogously for C3). The contribution from energies of O^m^) 
are usually split in three terms (actually this division is usually made irrespectively of the 
energy which is being integrated out): point-like, Zemach and polarizabilities corrections. 
They will be discussed further later. From the point of view of chiral counting the three 
terms are of the same order. Therefore, at the order of interest we can divide c^n R in the 
following way 

Ca^NR = ^A,R + ^Appoint-like + SCa, Zemach + ^'pol. ■ (35) 
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Indeed, a similar splitting is usually done for c\ l NR : 

C 3,NR = C 3,R + $ ^point-like + ^ c 3,Zemach + ^ c 3,pol. (36) 

Let us stress at this point that we are only interested in the logarithms. Therefore, we 
do not need to take care of the finite pieces. This will significantly simplify the calculation. 

We obtain the following result for the point-like contribution 

rvk 3 + 2c F - c| 2 m\ 

0C 4,point-like — ^ a 111 ~^2~ ■ K 6 1 ) 

Usually, in the literature, the computation of this type of contributions is made considering 
the proton point-like and relativistic, i.e. using standard QED-computations even at scales 
of 0(m p ). The fact that the proton has some anomalous magnetic moment that is due to 
hadronic effects or, in other words, the fact that proton has structure makes such theory 
no-renormalizable. This makes the result of the computation divergent. Within this phi- 
losophy, the result to the hyperfine splitting due to point-like contributions in Ref. [3] was 
proportional to 2 

3 + 2c F - c% 9 , mf. 3 „ m? 
£ L a 2 i n _|_ _ 

4 m 2 p 4 

where A is the cutoff of this computation. This computation would make sense if the scales 
on which the structure of the nucleon appears were much larger than the mass of the nucleon 
(then A could run up to this scale). Unfortunately, this is not the case and the structure of 
the nucleon appears at scales of 0(m p ) or even before. Therefore, to compute loops at the 
scale of 0{m p ) assuming the proton point-like produces problems (divergences) as we have 
seen. The procedure we use in this paper to deal with this issue is to work with effective 
field theories where the nucleon is considered to be non-relativistic. In other words, only at 
scales much smaller than the mass of the nucleon it is a good approximation to consider the 
nucleon point-like. Other usual method is to use some parameterization of the form factors 
fitted to the experimental data (see for instance [18, 19]). This regulates the ultraviolet 
divergences providing predictions for the hadronic correction to the hyperfine splitting. This 
is a very reasonable attitude in the cases where we are mainly interested in getting a number 
for the hadronic correction to the hyperfine splitting. Nevertheless, this is not the procedure 
we will follow in this paper, since our aim here is to gain as much understanding as possible 
on the structure of the proton from QCD and chiral symmetry. We want to understand 
how much of the coefficient can be understood from logarithms and energies of 0(m„-), for 
which a chiral Lagrangian can be used. This is something that could not be done with the 
standard form factors used to fit the experimental data, since they do not incorporate the 
correct momentum dependence at low energies due to chiral symmetry. 

The Zemach correction is due to what is called the elastic contribution, Eq. (27) of Ref. 
[16] for A 1 and analogously for A 2 (nevertheless A 2 does not appear to give contribution). 
It reads 

5cf Zemach = {^afm p \ J Jl^^GfG® . (39) 

2 Yet relativistic-type computations can be very useful sometimes, like in identifying some h\(mijmp) — > 
hi{mijv) logarithms and mijm p corrections in an efficient way, see [17]. 



(38) 
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Equation (39) can be obtained directly from Eq. (34) or working directly with the non- 
relativistic expressions and introducing the form factors. Either way, it is comforting to find 
the Zemach expression [1]. 

The upper index in and Gm has to do with the chiral counting. G$ = 1. It is 
illustrative to split the contribution to G$ from u and d, and the A: G$ = G^} d + G^ A 
and analogously for the Zemach contribution 

S^Zemach = ^,Zemach-u,d + ^,Zemach-A ■ (40) 

Another strong simplification comes from the fact that we are just searching for loga- 
rithms. Therefore, we are only interested in the behavior of the form factors for m p ^> k ^> 
m n and not analytical in k 2 . In particular, for the logarithms, we are only interested in the 
linear behavior in |k|. From Ref. [20, 21] (see also [22]), we obtain 



r (2) ^ ™ p 7T 2 

Gm ' a " (4^pl2 |k| 



(42) 



o 2 no, 



K,Zemach-u,d - (4^)2 a 3* 9a^~^, (43) 

- v l ^» 2^22 1 / a a \ 

K,Zemach-A - (4^)2 Q 27 ^ 9-kNA ln ^ ■ ( 44 ) 

It is remarkable that the above results are 7r-enhanced. 

Just for completeness, we also give the expression for S^ Zemach : 

^tzemacH = ^™fm\m h J ^^l^f Gf . (45) 

This term appears to be finite to the order of interest (it produces no logarithms) and it 
agrees with the result obtained by Pachucki [23] at leading order. 

The Zemach corrections (both the spin-dependent and the spin-independent) correctly 
incorporate the whole dependence on the lepton mass. Therefore, the same expressions are 
valid for the hydrogen and the muonic Hydrogen. 

Let us now consider the polarizability contributions. In the 577(2) case (and including 
the A) they should come from Eqs. (30,32,33,34,35,36) of Ref. [16]. In principle, in our 
case, one should consider more diagrams besides those plotted in Ref. [16], like the ones due 
to the Wess-Zumino anomaly action 3 . Nevertheless, as already stated in the introduction of 
the pionic Lagrangian, it turns out that they do not contribute in our case. Finally, we split 
the polarizability contribution as follows (for the SU(2) case): 

= 6 < P ol-A + K P ol.-,N + KU^A ■ (46) 



3 We thank T. Hemmert for stressing this possibility to us. 



10 



It is again a great simplification the fact that we are only searching for logarithms. 

From Eqs. (32,35) of Ref. [16] we obtain (we also checked this result by doing the 
computation directly in the non-relativistic limit) 



<o,-A = %« 2 ln^, (47) 



"h _ U 1,F 2 



where bf = G 1 according to the definition in Ref. [16]. The consequences this result has 
in a large N c analysis are remarkable enough. In the large N c , 6f = 3/ (2y/2) ny according 
to the Ji and Osborne definitions, where fiy stands for the isovector magnetic moment. On 
the other hand, by using the results of Ref. [24] \x v j \x n = —1, one obtains (for practical 
purposes) fj, p — /iy/2 in the large N c . It follows that in this limit the role of the delta is to 
cancel all the \i v contribution in Eq. (37) ((3 + 2cp — c 2 F )/A = 1 — /ip/4), which effectively 
becomes the result of a point-like particle. 

From Eqs. (30) and (34) of Ref. [16] we obtain 

8c? l < - ^ a 2 — -C In ^ (48) 

where C is defined in the Appendix. 

From Eqs. (33) and (36) of Ref. [16] we obtain 

c vi m l 2 " 2 64^, A 2 tAn . 

=(W 9 '-V27 Ch1 ^- (49) 

It is worth noting that Eqs. (48) and (49) cancel each other in the large N c limit, since g n NA — 
?>/(2\f2)g A in this case with the definitions of Ref. [16]. Moreover, they are suppressed by 
1/n factors and the smallness of the numerical coefficient compared with the Zemach term. 

Let us note that Eqs. (47), (48) and (49) may bring some light on why the polarization 
term is much smaller than the Zemach term in a model-independent way since we have an 
almost analytical result. 

Our results can be summarized in Eqs. (37), (43), (44), (47), (48) and (49). The above 
computation has been performed in SU(2), it would be interesting to repeat the analysis for 
SU(3). Indeed, we can compute the Zemach correction due to the strange quark (if we do 
not consider the spin 3/2 baryons) by using the results of Ref. [25] for Gm,s- 



g£) <* = (wi2 |k|[5Z}2 " 6DF + 9n (50) 



2 



We then obtain 

d D ~ l k 1 



K -(4^F a 1^ (5D2 - 6DF + 9F2)ln! ^- (51) 

In order to obtain a complete result, one should add the strange-related contribution to 
the Zemach correction due to the baryon spin 3/2 multiplet and to obtain the whole strange 
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contribution to the polarizability. This would require to have A\ and A 2 for SU(3) and 
including the baryon spin 3/2 multiplet, which are unfortunately unknown. In any case, one 
may wonder whether, for the polarizability corrections, the large N c cancellation would also 
hold in this case as well as the 1/ir and numerical factor suppression, so that it would be a 
very tiny contribution as in the SU(2) case. 

4.1 Matching to pNRQED: Energy correction 

With the above results one can obtain the leading hadronic contribution to the hyperfine 
splitting. It reads 

£hf = S^-{^ hp af . (52) 

By fixing the scale v = m p we obtain the following number for the total sum in the SU(2) 
case: 

^logarithms K) = "0.031 MHz . (53) 

The absolute value of this number would increase for a larger value of v and decrease for 
a smaller value. The main contribution to Eq. (53) comes from the Zemach and point-like 
corrections: 

^HF,Zcmach-u,d(m p ) = -0.022 MHz , (54) 
£HF,Zemach-A(m p ) = -0.004 MHz , (55) 
^HF,point-likcK) = -0.003 MHz. (56) 

Equation (53) accounts for approximately 2/3 of the difference between theory (pure 
QED) [18] and experiment [26]: 

E nF (QED) - E RF (exp) = -0.046 MHz. (57) 

What is left gives the expected size of the counterterm. Experimentally what we have is 
C4,nr = —48a 2 and c 4i _R(m p ) ~ —16a 2 . This last figure gives the expected size of the 
counterterm of the Lagrangian. A more detailed analysis would require to work in an specific 
scheme (for instance MS) to fix the finite pieces. We expect to come back to this issue in the 
future. A point to stress is that this number is universal, i.e. the same for the electron and 
for the muon up to corrections suppressed by the ratio of the lepton mass versus the proton 
mass: 

<*K) * c^(m,) . (58) 

This observation could be used in an eventual measurement of the hyperfine splitting of the 
muonic hydrogen. 

The introduction of the partial SU(3) computation would worsen the above prediction 

by 

^HF,Zemach-kaon(^p) = 0.003 MHz (59) 

bringing the total sum down to —0.027 MHz and c^R{m p ) to —20a 2 . 
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5 Conclusions 



We have performed a first exploratory study on the application of effective field theo- 
ries emanated from chiral Lagrangians on atomic physics 4 . We have computed the c^ l NR 
matching coefficient of the NRQED Lagrangian for the e-p and ji-p sector with 0(a 2 x 
(lnm g , In A, lnm/J) accuracy. The Hyperfine splitting of the hydrogen and muonic hydro- 
gen has been computed with 0(mf.a 5 jm 1 x (lnm q , In A, lnm^)) accuracy. We note that our 
results include the complete expression for the leading chiral logarithms. 

The difference between the experimental value of the hydrogen hyperfine splitting and 
the pure QED computation reads 

E HF (QED) - E HF (exp) = -0.046 MHz, (60) 

whereas our theoretical prediction reads 

^logarithms K) = -0.031 MHz . (61) 

We are then able to obtain an estimate of c\ R (m p ) ~ —16a 2 which is valid for both e-p and 
fi-p systems. 

One could improve these results by performing the whole computation in the MS scheme 
or alike (in this case some of the expressions in this paper should be rewritten in D- 
dimensions). It would follow that not only the logarithms but the finite pieces (in a specific 
scheme) of the matching coefficient would be obtained too. This would fix with a greater 
precision the value of c^ R (m p ) and, thus, the respective size of the effects due to the physics 
at scales of 0(m p ) and at scales of 0(m n ). This is important since the experimental number 
is precise enough to give an accurate number for c\ R (m p ), which could be used to test models 
at 0(m p ) scales. To perform the full computation in SU(3) would be also highly desirable. 
A partial SU(3) result brings Eq. (61) down to -0.027 MHz. 

Our results may help to better understand the fact that the Zemach correction is much 
larger than the polarizability contribution, since we have (almost) analytical expressions for 
these contributions. The polarizability term (except Eq. (47)) vanishes in the large N c and it 
is 1 /V and numerical-factor suppressed with respect the Zemach terms. On the other hand, 
Eq. (47), in the large N c limit, cancels all the fi p contribution in Eq. (37), which effectively 
becomes the result of a point-like particle. 

Several lines of research are worth pursuing. One is trying to compute c^^r within 
HBET, since its numerical value could be obtained from measurements of the Lamb shift, 
and it is related with (and in a way defines) the proton radius. Another could be to consider 
more complicated atoms within this effective field theory formalism (see, for instance [9] for 
the Helium). 
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Appendix A Constants 



T— 1 

^0 


= 92.5 MeV, 


9a 


1 c 
= l.zo , 


m w 


= 14U Me V , 


m p 


= m n = 938 MeV , 


A 


= 294 MeV, 




= 1.05, 


b( 


= 3.86, 


F 


= 1/2, 


D 


= 3/4, 


m p 


= 770 MeV. 



(62) 



bf and g^NA have been obtained from the decays of the delta in the non-relativistic limit 
(consistent with the accuracy of our calculation). 

The values of F and D are consistent with the large N c limit. Finally 



y \ v ' V x ~ x + x V 1 



- 3 < J - 2x) y 2 i i- x l- f) ) Sinlrl (\/(r^) v) j dvdx 

= -0.165037. (63) 
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